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Abstract 

Having in view some applications in nanophysics, in particular in 
nanophysics of materials, we develop new dynamical models of structured 
bodies with affine internal degrees of freedom. In particular, we construct 
some models where not only kinematics but also dynamics of systems of 
affine bodies is affinely invariant. Quantization schemes are developed. 
This is necessary in the range of physical phenomena we are interested in. 

Introduction 

The idea of micro-structure is a rather old one and goes back to brothers Cosserat 
who formulated the theory of continuum consisting of infinitesimal gyroscopes 
[1 ]. Eringen modified such a model by introducing homogeneous deformations 
as additional microstructural modes [2] ■ These theories were rather phenomeno- 
logical and often motivated by some kind of mathematical (more precisely, 
differential-geometric) aesthetics. There was also mechanical motivation based 
on theory of granular media and continuous limit of the dynamics of molecular 
crystals. 

Recently the interest in mechanics of structured media becomes more and 
more intensive in connection with nano-structures, various supramolecular stru- 
ctures, and defect theory. There are also some special problems like the dynam- 
ics of suspensions, gas bubbles in fluids, and some very peculiar models like 
kinetic media [3] HI [H [7] • In all these problems, in defect theory and in dy- 
namics of fullerens, affine model of structural elements is well- motivated both 
from the physical and geometrical point of view. There are some interesting 
problems, completely new ones in comparison with traditional phenomenolog- 
ical mechanics of structured continua. First of all, in dynamics of strongly 
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interacting systems of objects in nano-scale, on the classical level, the tradi- 
tional scheme of constrained motion based on the d'Alembert principle does 
not seem reliable any longer. Instead, one should base the dynamics on some 
effective collective models motivated by appropriate symmetry demands. We 
believe, there are physical reasons to expect that the dynamics of internal and 
collective affine modes should be also affinely invariant. In traditional theories 
derived from the d'Alembert principle, the dynamics of affine modes is invariant 
only under the Euclidean group. Below we suggest some dynamical models of a 
system of bodies with affine degrees of freedom, with Hamiltonians invariant un- 
der the affine group. One can smoothly include also some terms invariant only 
under the Euclidean group, there is however a natural temptation to consider 
them as a merely perturbation to the background affine dynamics. Some argu- 
ments from the solid state physics, like the concept of effective mass support the 
idea that in complicated problems of condensed matter theory the "true" metric 
tensor is not necessarily a fundamental geometric background of equations of 
motion, Lagrangians and Hamiltonians (having in mind microphysical models, 
we assume Hamiltonian, variational dynamical models, without dissipation on 
the fundamental level). 

There is also another important novelty in comparison with traditional 
macroscopic models. Namely, in the nano-scale the quantum background of 
the dynamics must be seriously taken into account. Because of this we devel- 
oped some quantization scheme. But at the same time one uses the concepts of 
macroscopic origin, like, e.g., deformation tensors, deformation invariants, etc. 
It seems rather amazing that one formulates questions like: what is a complex 
quantum probability amplitude for the deformation tensors and deformation in- 
variants to be found in some fixed range? What are the corresponding quantum 
transition probabilities? There are even more serious problems. Namely, in this 
very peculiar range of phenomena one has to do with a very complicated convo- 
lution of classical and quantum problems. In this, rather unexpected way, the 
old problems from the realm of foundations of quanta [S], like decoherence, wave 
functions reduction, possibility of nonlinear quantum description, etc. revive as 
ones motivated by quite practical, structural physics. One must honestly say, 
there is more secrets than well-established facts and answers here. 

1 Affine bodies and their systems. General con- 
cepts and basic quantities 

There is no place here for the very detailed geometric description. A rather 
exhaustive treatment of the differential-geometric background of our ideas was 
presented in some earlier papers (cf. [10] and references there) and a more 
detailed treatment will soon appear. Here we base mainly on the analytic de- 
scription where the physical and material spaces are simply identified with R n 
by some choice of orthonormal Cartesian coordinates. As far as possible we work 
in a non-specified dimension n and only at some final stages we specify n to its 
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physical value 3 or to 2 and 1, which obviously are also physically interpretable. 
This "false" generality is mathematically convenient and better reveals some 
structural features of the model, hidden behind the particular value n = 3. 
Moreover, it suggests even some analytical solving procedures. 

So, the body we consider consists of elements performing the translational 
motion in the physical space R™ (physically n = 3,2, 1) and also some internal 
motion. The latter is simply a relative motion of microconstituents; the trans- 
lational R"-degrees of freedom are attributed to the centres of mass. But of 
course one can also admit situations when additional degrees of freedom are es- 
sentially internal ones, like, e.g., in spin media. And in any case the distinction 
between internal and relative motion may be (although need not be) historical, 
based on some conventions or on our laboratory abilities. From now on we do 
not discuss this problem and for brevity all non-translational degrees of freedom 
will be referred to as internal ones. 

In mechanics of structureless continua the configuration space may be iden- 
tified with Diff(n,R), i.e., the group of (sufficiently smooth) diffeomorphisms of 
R™ onto itself. Let us remember, according to our analytical conventions both 
the physical and material space are simply identified with R™. Any diffeomor- 
phism ip G Diff(n, R) establishes an interrelation between Lagrangian (material) 
and Eulerian (physical, current) coordinates, respectively a k and x 1 , 



In this way the configuration space of continuous medium is identified with the 
infinite-dimensional group Diff(n,R). In the Arnold description of incompress- 
ible ideal fluid [H] this group is constrained to SDiff(n,R), i.e., the subgroup 
of volume-preserving diffeomorphisms. Euler equations are interpreted in terms 
of right-invariant geodetic Hamiltonian systems on this group. This approach 
turned out to be at least hcuristically effective in hydrodynamical problems. 

Of course, Diff(n, R) would be completely non-effective and just meaningless 
as the configuration space of internal degrees of freedom of structured bodies. 
Nevertheless, there exist both physical and analytical reasons to concentrate on 
models with degrees of freedom based on appropriately chosen groups and their 
homogeneous spaces. From the physical and computational point of view it is 
clear that one must use some geometrically well-motivated finite-dimensional 
subgroups G C Diff(n,R), just Lie groups acting in R ra . The most traditional 
pattern is G = SO(n,R), i.e., rigid body model of internal degrees of freedom. 
It is not excluded, especially in quantized theory that the non-connected con- 
figuration space G = 0(n,R) with mirror-reflected configurations may be also 
acceptable. Other natural model with more degrees of freedom is that of affincly- 
rigid, i.e., homogeneously deformable, body; G = GL+(n,R) or perhaps, as 
above, the total non-connected linear group G = GL(n,R), with reflections ad- 
mitted. One can also think about L(n, R), the algebra of all n x n matrices as an 
admissible configuration space. Incidentally, there is a subtle difference between 
admitting mirror- reflected configurations of metrically- and affinely-rigid bodies. 
Namely, GL + (n,R) and its mirror-reflected coset (not a subgroup!) GL~(n,R) 
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in GL(n,R) infinitcsimally approach each other in L(n,R), being separated only 
by the (n 2 — l)-dimcnsional subset of singular matrices in L(n, R). On the other 
hand, SO(n,R) and the set of improper rotations 0(n,R)\ SO(n,R) (the com- 
plement of SO(n,R) in 0(n,R)) are so-to-speak finitely separated in L(n, R). 

There are also natural models placed between gyroscopic and affine degrees 
of freedom; one deals then with a constrained affinely-rigid body. Let us mention 
incompressible (isochoric) body when G = SL(n, R), the special linear group. 
Just as previously, one can also think about admitting mirror reflections. Then 
G = UL(n,R) = {ip G GL(n,R) : |det(p| = 1}, the unimodular group. It is a 
union of SL(?i, R) and of its coset in GL(ro,R) consisting of matrices with the 
minus-one-determinant. Obviously, just as in the case of metrically-rigid body, 
the manifolds of proper and improper (orientation-changing) isochoric mappings 
are finitely-separated in L(n,R). In a sense an opposite model (less realistic) 
is that of shape-preserving affinely-rigid body when, roughly speaking, internal 
configurations are built of dilatations and rotations, G = R + SO(n, R); R + 
denoting the multiplicative group of positive real numbers. Again, admitting 
orientation-preserving mappings we have G = R + 0(n,R). And here also the 
two connected components approach each other infinitely close. The closure of 
G contains the null matrix, therefore, it is not a subset of GL(n, R). 

All these subgroups of Diff(n,R) are contained in GL(n,R); they consist of 
linear mappings. These mappings describe configurations in such a way that 
the material point with Lagrangian coordinates a k occupies the spatial position 
with Eulcr coordinates y l , 

V 1 =x l + <p l k a k , (2) 

where x 1 are spatial coordinates of the centre of mass. The very geometry 
suggests also some groups of non-affine transformations, e.g., the projective 
group GPr(n, R), conformal group Co(n, R). One can realize physical situations 
where the relevant collective modes of internal motion are described just by these 
groups. It is also non-excluded that the complex group GL(n, C) or its unitary 
subgroup U(n) C GL(n, C) 12J may be useful. In any case, it happens quite 
often in physics that the complexification leads to quite amazing, unexpected 
new results. 

Just as in our earlier papers we concentrate here on the affine model, when 
G = GL(n,R). Incidentally, the projective model may be in some sense re- 
duced to it, because GPr(n,R) is isomorphic with SL(n + 1,R). In any case, 
affine modes seem to be dominant for small structure entities like molecules, 
microdefects, some supramolecular clusters, fullerens, etc. 

The configuration space of a single structural elements is identified with 

Q = M" x G = Q tr x Q mt , (3) 

in particular, for elements with affine modes of deformation: 

Q = R" x GL(n,R). (4) 

Usually, especially in classical (non-quantized) problems, GL(n,R) is replaced 
by GL + (n,R). The labels "tr" and "int" refer obviously to transnational and 
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internal degrees of freedom. The total body (medium) consists of N elements. 
Its configuration space is obviously given by the Cartesian product 



in our treatment GL(n,R) or GL + (n,R) substituted for G. Therefore, the 
configuration is an array: 

q = (x!, . . . ,x N ;ip!, . . . ,tp N ) , (5) 



where x A G R™, ip A £ GL(n,R), and A = 1,N. 

For a single structure element the summation of usual kinetic energies of its 
constituents gives in virtue of ([2|), the usual d'Alembert form 

T d'A = T d<A + T £A = {vyT) + 1 ^ (g) 

where v € R n , £ £ L(n, R) denote respectively the translational and internal 
velocities: 

and M, J are constant inertial characteristics. More precisely, M is the total 
mass of the element and the symmetric positively definite matrix J is its mod- 
ified inertial tensor, i.e., second-order moment of the mass distribution with 
respect to co-moving (Lagrange) coordinates, 
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Summation is performed over constituents ("atoms") of the element ("mole- 
cule"); fi p is the mass of the p-th constituent. Sometimes it is convenient to use 
the symbol of integration with respect to the mass distribution measure \x: 



J kl = J a k a L d^{a). (9) 

Remark: the kinetic energy ([6]) is spatially isotropic, i.e., invariant under 
the transformations La below (fT4|) with A restricted to the orthogonal group 
0(n,IR) (spatial rotations). So are its both terms separately. The material 
rotations Ra in (]14[) preserve T tr trivially, but in general Ti nt is non-invariant 
under the right-acting 0(n,R). However, it is invariant under the right actions 
of 0(n, J), the subgroup of GL(n, R) preserving J. This reduces to the 0(n, R)- 
invariance, when J is isotropic, i.e., J = /Id„; / is a positive constant of internal 
inertia and Id rl is the n x n identity matrix. Then 

7?„; A ^Tr(^). (10) 
The total kinetic energy of the body is given by 

N N N 

T dA = J2 T a A = \Y. M ^ M) + \ £ Tr {U J Ail) ■ (11) 

A=l A=l A=l 
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Assuming that the body consists of identical structure elements we have that 
M A = M, J A = J, A = TJ?. 

Let us again concentrate on a single element. Its Green and Cauchy defor- 
mation tensors are respectively denoted as 

G[<p] = <p T <p = G[^] T , C[<p] = tp-^tp- 1 = CM T , (12) 

similarly, for their contravariant inverses we write 

G[<p] = <p- 1 'p- 1T , CM = w T . (13) 

Spatial and material transformations are respectively given by left and right 
regular translations: 

ip ^ L A (ip) = A<p, tp h-> R A (<p) = ip A (14) 

for any fixed A £ GL(n,K). When A E 0(n, K), then obviously 

G[A<p]=G[<p], C[ipA)=C[pl (15) 

and for the general A G GL(n, R) 

G[tpA] = A T G[(p]A, C[Ai P ]=A- 1T C{ip]A- 1 . (16) 

There is no concise formula for G[-A<^], C[y>A] if A is not orthogonal (does not 
belong to 0(n,M)). 

Deformation invariants are scalar functions / : GL(n,R) — > M invariant 
under (jTJJ) for orthogonal translations 

f(A<pB) = f(<p) (17) 

for any A,Be 0(n,R). There are n basic invariants through which all other 
ones may be expressed. Various choices are possible, e.g., the following fre- 
quently used 

Ka[<p] = Tr (G[cp] a ) = Tr (C[<p]~ a ) , a = T/n, (18) 
eigenvalues X a [f] of G[ip], 

det (G[tp] - X[ip]I n ) = 0, (19) 
or coefficients I p [ip] of the eigenequation 

n 

det (G[<p] - XIn) = £(-l)*/„_fcMA*; (20) 

fe=0 

obviously, Iq = 1 is standard. Geometrically speaking, deformation invariants 
are functions on the manifold of double cosets 

Inv := 0(n,R)\GL(ri,R)/0(n,R). 
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Deformation invariants are used when constructing potential energy models for 
a single affine body. When dealing with the system of such bodies we need some 
basic scalars assigned to pairs of internal configurations. In analogy to Green 
and Cauchy deformation tensors for any pair ip,ip € GL(n,R) we define the 
quantities 

G[i>, ip] := C[V, <p] ■= ?- XT ^-\ (21) 

Obviously, 

G[il> i ft>]=G[fl>], C[^,V]=CM, 

i.e., the above mutual deformation tensors reduce then to the usual ones. 
But one can also define another mutual quantities, namely, 

For orthogonal matrices they reduce to the previous ones, 

^eO(n,R)=*rfoM=GfoM, T.[1>,<p\ = C[il>M- (23) 

Obviously, r[^>, if], Tl[ip,ip] are exactly group-theoretical counterparts of the 
displacement vector in translational degrees of freedom. Indeed, interpreting 
R™ as an Abelian group under addition of vectors, we immediately notice that 
the prescription (|22p in the non- Abelian multiplicative matrix group GL(n, R) 
has exactly the same group meaning as u = y — w in R™. 
It is clear that for any A £ 0(n, R) 

G[Ail>, A<p] = Gfof, f], CtyA, V A] = C[il>, <p], (24) 

i.e., they are respectively invariant under spatial and material isometries. For 
the general A £ GL(n, R) we have 

G[ipA, f A] — A t G[i{j, f]A, C[Aij),Aip\=A- 1T C[^,y\A- 1 . (25) 

And just as previously there is no concise expression for G[Aip, Atf], C[ipA, ifA] 
if A is non-orthogonal. 

Transformation rules for T, S have another form. Namely, for any A G 
GL(n,R) we have 

T[Aif>, Asp] = T[i>, if], V[A1>, Asp] = AY,[^, f]A~\ (26) 

TtyA, fA] = A- l T[^, ip]A, m>A, if A] = E[^, if]. (27) 

Therefore, T is invariant under spatial affine transformations and suffers the 
inverse adjoint rule under material affine transformations. And conversely, E 
transforms according to the adjoint rule under spatial affine mappings and is 
affincly invariant under material transformations. 

The quantities G[ip,if], C[ip,if], T[ip,if], Y,[ip,if], give rise to scalars which 
may be used as arguments of the potential energy terms. Typical scalars of this 
type are, in analogy to (|18p . given by 

/C a ^,^=Tr(G^,^] a )=Tr(C^,^]- a ), a = T~n. (28) 
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Just as in the case of deformation invariants, these scalars are invariant under 
spatial and materia rotations (left and right regular translations of ip, ip by 
orthogonal matrices) : 

K a [Ai,B,A<pB]=K a [ip,<p}, A,BeO(n,R). (29) 

In analogy to (fT9|) . (|20| one can also use solutions of the eigenequation for 
G[ip, if] (or C[ip, if]), or coefficients in the eigenequation as basic invariants. 
Another kind of invariants is built of T, E-objects, e.g., 

M a [i/>M =Tr(r[^] Q ) = Tr(£[^f) (30) 

like in (fl%)) . or, according to the A a -, I a -schemes like in (TT9"1) . ([2l)|) . These objects 
are invariant under all affine spatial and material transformations, i.e., 

M a [AipB,AipB}=M a [ip,ip}, A,fl6 0(n,R) (31) 

for any A,Be GL(n, R). These scalars measures of the "distance" between 
internal configurations are affinely invariant. Unlike this, the measures IC a are 
only orthogonally invariant, so they are usual Euclidean distances. 

Let us remind that in many problems it is convenient to use deformation 
measures which vanish in the non-deformed state, e.g., Lagrange and Euler 
deformation tensors: 

E[<p] = ~ (G[ip] - I) , e[ip] = \ (I - C[ip]) . (32) 

By analogy, it may be convenient to replace the mutual tensors G[ip,cp], C[i/),tp], 
r[ip,(p], T{ip,<p] by 

%*>] = 5(GM-i), e[^M = \{I-CVp\), (33) 

~/ty,tp}=Tty,<p\-I, tr[1>,tp]=m ) ,<p]-I- (34) 

Another possibility is to use the exponential representation of matrices G, C, 
r, S. It may be also convenient to use invariants built of them according to 
the schemes like (|T5|). (1501) or (TT9")) . (J2UJ), etc. Obviously, such invariants are 
functionally dependent on the previous ones. 

Affine velocity, i.e., Eringen's "gyration" [5j respectively in the spatial and 
co-moving representations is given by 

Sl=%-\ fi = *T* (35) 
dV Y dt K J 

Spatial and material transformations (TTi)) act on the above quantities as follows: 
L A : Q^AQA' 1 , f2 >->• Q, (36) 

r a : n >-> o, n >-> a^vla. (37) 
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When the motion is metrically rigid, i.e., permanently ip 6 0(n, K), then 

n = -o T , h = ~h T (38) 

and these skew-symmetric objets reduce to the usual angular velocity, respec- 
tively in the spatial and co-moving representations. In some formulas it is 
convenient to use also the co-moving representation of translational velocity. It 
is given by 

v^p^v. (39) 

Gyroscopic constraints may be described in anholonomic terms simply by 
stating that Q, is skew-symmetric (and so is VL then), 

o + n T = o. (40) 

By analogy, constraints of the purely deformative rotationless motion may be 
defined by the demand of permanently symmetric f2, 

n - n T = o. (4i) 

Let us observe that the materially rotationless constraints 

- h T = (42) 

are non-equivalent to the above spatially rotationless ones. 

There is an interesting novelty now. Namely, rotationless constraints are 
essentially non-holonomic. The point is that the subspace of symmetric matrices 
is not a commutator Lie subalgebra. 

Canonical momenta, i.e., dual objects of velocities (linear functions of them), 
p, 7r are elements of R n , L(n,R) respectively, and their pairing with velocities 
is given by 

(p,v)=p T v = Tr(pv T ), (7r,£>=Tr(7rO. (43) 
The co-moving representation of p is given by 

p = <p T p, (44) 

and obviously, with this convention 

(p, v) = p T v = p T v = (p,v). (45) 

It is convenient to introduce non-holonomic canonical momenta conjugate 
to f2, O, namely, 

S := (pir, E := mp. (46) 

Obviously, 

(£,fi) = Tr(£f2) =Tr(lifi) = (47) 
Transformation rules for S, S are identical with (l36l). (l37l): 



L A : S i-» A'SA^ 1 , S^E, (48) 
R A : E^E, E^A^EA (49) 
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The components of £, £ are respectively Hamiltonian generators of spatial and 
material affine transformations (|14p . Their doubled skew-symmetric parts 



S = £ - £ J = -S, V = £ - £ J = -y (50) 

are, respectively, Hamiltonian generators of spatial and material rotations of 
internal degrees of freedom. 

The relationship between Hamiltonian quantities like p, p, £, £ and kinemat- 
ical ones like u, v, £1, H may be established only on the basis of some particular 
model, when Lagrangian L is fixed. In potential models L = T — V with the po- 
tential V depending only on the configuration (x, <p) G M™x GL(n,R), Legendre 
transformation has the following form: 

Assuming the d'Alembert model of the kinetic energy we obtain that 

p = Mv, n = Jf, (52) 
and the corresponding expression for the kinetic Hamiltonian 

T = T tr + 7I nt = ^Tr (pp T ) + ^Tr (tt t J^tt) . (53) 

For Hamiltonian systems H = T+V (V denoting the potential energy) equations 
of motion may be effectively analyzed in terms of Poisson brackets and canonical 
Hamilton equations, 

^ = {F,H}, (54) 

where F runs over some maximal system of functionally independent phase- 
space functions. 



2 Dynamical models. Affine invariance prob- 
lems. Realistic questions, academic questions, 
and pure fantasy 

For the system of affine bodies Lagrangian has the form: 

L = T-V, (55) 

where the kinetic energy is obtained by summation of individual kinetic energies 
like in (jlT]) . 

N N N 

T=Y, T A= T tI + T int = £ (T tI ) A + i^nt) A • (56) 

A=l A=l A=l 
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The potential energy in typical situations consists of two main terms, the one- 
and two-body potentials, 

V = V^+V^ 2 K (57) 

It is known that in realistic problems it is usually less then 10% of energy that 
could be assigned to three-body and higher multibody interactions. is the 
sum of terms depending on individual elements, 

AT 

V« (. . . ; x A , <p A ; ■ ■ •) = Yl V(1) B (xb,Vb) ■ (58) 

The over-simplified models where b splits into the sum of translational and 
internal parts, 

V« B (x B , <p B ) = V^B (x B ) + V^B {<p B ) , (59) 

are not very realistic, nevertheless, they provide some so-to-speak zeroth-order 
approximation. When the elements are identical, all b have the same func- 
tional form. 

The binary term has the usual form, 

1 N 

V (2) (. . . ;x a ,Va; ■ ■ • ) = 7r X! v( ~ 2) kl (x K ,ip K ;x Ll ip L ) . (60) 

K,L=1 

And again the simplest, although rather academic models are those with the 
separated dependence of on translational and internal variables, 

V^ 2) kl{xk,vk]x l ,lp l ) = V t (2) kl(xk,x l ) + V^kl('Pk,'Pl) ■ (61) 

Mutual interactions should be translationally invariant, i.e., V^ 2 ' kl depend on 
xk, xl through xkxl = xl — xk- Isotropy of the physical space implies that 
the radius-vectors ij, — xk enter kl only through their lengths \\xl — xk\\- 
There is some more discussion concerning the dependence of V' - 2 ' kl on internal 

(2) 

degrees of freedom. Isotropy of the physical space implies that V int kl should 
depend on ipK, fL only through the mutual tensors G[<£k, <Pl], F[<Pk, <Pl], thus, 

V {2) kl (x k ,tp k ;x L ,(pL) = !kl {\\x l - x K \\ ,G[ipk,<Pl] ,T[pk,¥l\) (62) 

and obviously for the body consisting of identical elements there is no depen- 
dence on K, L; /kl = f for some fixed /. And if the dynamics is to be invariant 
also under simultaneous material rotations, then at the same time, kl must 
depend on internal configurations only through C[(fK, <Pl], ^[<Pk, <Pl]- But this 
means that kl is algebraically built of the mutual invariants, e.g., chosen 
as lC a [ip K ,ipL], M a [ipK,¥>L], 

V (2) kl (xk, tpK;x L ,ifL) = Ikl [\\xl - X K \\ ,^[V>K,V>l] ,M [pk,Vl]) ■ (63) 

In the last formula, JC, M arc abbreviations for the systems K, a , M a , a = 1, n. 
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In our model, geometry of degrees of freedom and kinematics is ruled by 
the affine group. On the other hand, the dynamics is not invariant either un- 
der spatial or material affine transformations (fH)) . The spatial metric tensor 
and the inertial moment J break the affine symmetry and restrict it to the Eu- 
clidean one in the physical space and to 0(n, J) in the material space. What 
concerns potential energy of mutual interactions it is clear that the vector 
norm \\xl — %k\\ an d transposition-dependent invariants K, [-Pk, <Pl] also restrict 
the spatial symmetry to 0(n, R). But it is well-known that particularly inter- 
esting models and successful analytical procedures appear when the group of 
dynamical symmetries (symmetries of Lagrangian) coincides with the kinemat- 
ical group, or at least, when it is as large a subgroup as possible. The questions 
arise as to the formal possibility and physical usefulness of affinely-invariant 
models. For a single affinely-rigid body such models are in a sense possible 
O [10]. Their physical usefulness is not yet decided, although there are some 
arguments supporting it. Namely, it is quite possible that in complex media 
with a complicated net of internal interactions a single element is more sen- 
sitive to its material surrounding than to the "true" metric tensor (produced, 
according to General Relativity by the gravitational field as its "vacuum" non- 
excited state). The more so such a mechanism works in defect theory. Let us 
also mention the concept of effective mass [T31 H3] in crystals, where the kinetic 
energy of electrons is not based on the "true" metric, but on the effective tensor 
produced by the material surroundings. There are nice mathematical models 
of the kinetic energy of a single affine body with the kinetic energy based on 
the Cauchy tensor used as a metric. There are also some physical arguments 
supporting such a hypothesis [TO] . 

The material affine invariance may seem perhaps more natural because there 
exist models of continua based on very rich material symmetry. As mentioned, 
this is Arnold description of the ideal incompressible fluid. It is based on infinite- 
dimensional group of volume-preserving diffeomorphism. They act on the right, 
i.e., as material transformations. In any case, finite-dimensional geodetic models 
of small grains or suspensions with kinetic energies materially invariant under 
SL(n, R) may be considered as on over-simplified, drastically discretized version 
of the Arnold model. 

Obviously, apparently the Euclidean invariance of the kinetic energy of single 
elements and of mutual potential elements seems to be firmly established. But 
nevertheless it may be a superficial illusion and affine invariance should be 
at least admitted to consideration. When one deals with a highly condensed 
matter and with very small structure elements, e.g., in the nanoscale, then the 
complicated structure of interactions may result in quite unexpected results. 
Namely that some hypothetic phcnomenological models based only on some 
symmetry guiding hints may be so realistic as (or even more realistic) than ones 
based on apparently careful structural " derivation" . It is so because " derivation" 
in strongly interacting structured media always neglects a lot of factors and the 
collective effective phenomena are better described by phenomenological models 
derived on the basis of well-established invariance assumptions. It is so, e.g., in 
such complicated structures like atomic nuclei. In any case, symmetry principles 
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arc then rather reliable guiding hints. 

For a single affine body non-trivial potentials of internal degrees of freedom 
are never affinely invariant. Only constant functions on GL(n, R) may be so. 
The same is true for one-particle external potentials of the system. As we have 
just seen, for the purely mutual interactions, the binary potentials of internal 
degrees of freedom admit affine invariants as arguments. As we shall see, at 
least formally, the same is true for translational degrees of freedom. We shall go 
back to this problem later on and concentrate now on the kinetic energy terms. 

There exist at least academically interesting kinetic energies for single affine 
bodies, and obviously, for their systems as well (because, unlike potentials, the 
kinetic energy is additive). 

Of course, the usual d'Alembert model ([6]) of Tj nt is isotropic in the physical 
space and affinely invariant under material affine transformations. This is also 
explicitly visualized by its another equivalent representations, 

Tlr af = ^Tr (vv T ) = ^v T v = ^v T G [ip] v = ^Tr (vv T G [<p]) . (64) 

And the same holds for any of the structural element, i.e., for any Ta- Obviously, 
v, v, G must be then replaced by the corresponding va, va, Ga — G[(Pa]- 

The labels "tr", "is", "af in (|6~4"|) refer respectively to "translational part", 
"isometry invariant in the physical space" (thus, written on the left-hand-side: 
left-invariant), and "affinely invariant in the material space of a single element" 
(thus, written on the right-hand-side: right-invariant). 

From the purely academic point of view one can also think about T^~ ls , 
the model affinely invariant in the physical space and isometry invariant in the 
material space (respectively, left and right invariance in GL(n, R)). It will have 
the form: 

T t 1 - is = fv T v = f Tr (vv T ) = fv T C [<p] v = ^Tr (vv T C [<p]) . (65) 

The Cauchy deformation tensor is now used as the "metric" of the physical 
space. 

Let us observe, there is no possibility to obtain translational kinetic energy 
which would be affinely invariant in both the spatial and material sense, at least 
if we do not try some extremely exotic things. The reason is that the affine group 
is non-simple in a very special way and does not admit doubly-invariant and 
non-degenerate twice covariant tensor fields. Let us observe that the trick with 
the Cauchy tensor substituted instead of the "usual" metric tensor may be as 
well repeated for the internal part of ©. This would lead to the following 
expressions: 

T t f- J = ^Tr (C [if] HJe) = ^Tr (hjh T ) . (66) 

Here the upper-case labels at T; nt mean that it is invariant under all spatial 
affine transformations and under the material group 0(n, J) C GL(rt,R) of J- 
preserving material transformations (" material isometries" when interpreting J 
as kind of the "material metric tensor"). The expression (|66[) becomes T-^ t ~ ls , 
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i.e., materially isotropic in the usual sense when J is spherical, i.e., J = IlA n , 
where I is the scalar inertial parameter and Id„ is the unit n x n matrix. Then 
we have 

T t f" is = T -Tr (nn T ) = T -Tr (C [<p] ■ (67) 



Let us observe that the spatially affine models (|65|) , (|66|) are constant- 
coefficient quadratic forms only when expressed in terms of non-holonomic ve- 
locities f2, therefore, the metric tensors underlying them are curved, essentially 
Riemannian. 

One can wonder what would be a possibility symmetric to (|66p . i.e., affinely 
invariant in the material sense. In this sense the natural modification of Tj nt in 
([6]) would be 

T«- af = ilr (Min T ) = llr (£H [<p] , (68) 

where H is a fixed positively definite matrix, and H[ip] is its following ip^ 1 - 
transform: 

H[ip}=ip- 1 Hip- 1T . (69) 

So, from the point of view of the usual d'Alembert theory H [ip] is a strange 
(^-dependent although co-moving inertial tensor, whereas its spatial represen- 
tation is constant. In the usual d'Alembert affine dynamics the co-moving J 
is constant, whereas its spatial representation J[ip] is configuration-dependent, 
thus, also time-dependent: 

J[ip\=ipJip T . (70) 

The above affine models of Xi nt are very peculiar in the sense that the inertial 
terms in (|6"6"]) , (f6"8")) are factorized into tensor products. (The factors additional 
to H, J are the physical and material metric tensors. They are apparently ab- 
sent because we use Cartesian coordinates which analytically reduce the metric 
tensors to identity matrices.) The most general Tj nt invariant under spatial and 
material affine transformations are respectively given by 

^ af - \c B A D c n A B h c D , (7i) 

= l -WhW 3 n k u (72) 

where £, 1Z are constant and symmetric in their bi-indices (under exchanging 
mutually the first and second pairs of indices). 

The only model of Ti nt invariant simultaneously under spatial and material 
affine transformations (left and right regular translations in GL(n, R)) has the 
form: 

7f„ f r af = ^Tr (f7 2 ) + f (Tr nf = ^Tr (fi 2 ) + | (Tr n) ' , (73) 

where A, B are inertial constants. 

This affine-affine model of Tj nt is never positively-definite. One can show that 
this "failure" is non-embarrassing, moreover, it may be just profitable [9| flO] . 
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It was told above that the translational part T tr is never simultaneously left 
(spatial) and right (material) afhnely invariant. The highest available symme- 
tries of Ttr are T t ^~ af and T t a r f " is f[H5 ]> . This raises the question as to the 
internal counterparts T-^ a and T^ -153 . One can easily show they are given by 



where / is an additional inertial constant. In the special case of the metrically- 
rigid (gyroscopic) motion, the first terms in ([74]) . ([75)) coincide. In general this is 
not the case. Therefore, if gyroscopic constraints are imposed (no deformations), 
(j74[) and (Tf5|) coincide, and one obtains the spherical rigid body with the scalar 
inertial momentum (I — A). 

Unlike ([75)) . expressions (|74p. ([75]) may be positively definite and they are so 
in some open range of triples (/, A, B). And at the same time they continue to 
have the main advantages and the general structure of (|73[) . both on the level 
of theoretical analysis and practical calculations. 

Obviously, everything said above concerns directly a single affinely-defor- 
mable element, nevertheless, just as in the d'Alembert model, applies also im- 
mediately to the total system, because the modified kinetic energies are additive. 
One should only use explicitly the label K referring to structural elements. 

And now let us go back to the problem of potential energy. As mentioned, the 
external one-particle potential V*- 1 ^ cannot be affinely invariant (only constant 
functions may be so) . The formula (|63p for the doubly isotropic binary potential 
V KL seems to suggest something similar for the dynamics of mutual interactions. 
However, things are not so simple and one can try to find some modifications 
towards the affine invariance, just as it was done in the case of kinetic energy 
models. Some at least formally admissible suggestion may be easily formulated. 

Let us fix some pair of structural elements labelled by (K,L). Internal 
configurations (Pk,<Pl £ GL(n, K) give rise to the Cauchy tensors C^k), C[vl]- 
In our considerations above we were faced with the idea of using C[tp] as a kind 
of spatial " metric tensor" underlying affinely- invariant kinetic energies of single 
elements. Let us now introduce the objects 



It is symmetric in the labels K, L and positively definite. This motivates the 
temptation to use it as a "metric tensor" underlying some modified "distance" 
between xk and xl, namely, 




(74) 



(75) 



C[<p K> <p L } = -(C [<p K ] + C[<p L ]). 




(76) 
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Obviously, the transformation rule (fill)) implies that the above prescription is 
invariant under the spatial action of GL(n, R): 



V [Ax K , Aip K ; Ax L , Aip L ] = V [x K , ; x L , tp L ] (77) 

for any A S GL(n, R). This is a rather curious affinely-invariant "distance" . And 
now we can modify (|63p by introducing to it this new distance-like argument in 
addition to the usual one: 



\\xk - xl || = \J {xk ~ xl) T {xk ~ xl) = J Tr ([x K - xl) {xk - xl) T ^ ■ 
So, finally, instead of (f63|) we have 



(2) 

v kl ( x k,<Pk;x l ,<p l ) = (78) 
= !kl [\\xk - x L \\ ,T>[x K , Pk;xl,<Pl] ,JC[fK, ¥l] ,M [<Pk,<Pl]) , 

where in realistic situations all /kl coincide with some fixed /. It is seen 

(2) 

that V KL depends on its configuration arguments through the system of four 
scalar quantities. Two of them, namely, \\xk — and K, [<pk, Pl] are invari- 
ants of the rotation group 0(n, R). The remaining two, T> [xk, <Pk', xl, <Pl] and 
M. [ifK, Vl], are invariant under the total linear group GL(n,R). One can ex- 
pect that the dependence of on the latter two scalars is a highly symmetric, 
affmc background of mutual interactions between constituents of the body. Fur- 
ther on, this high affine symmetry is broken and reduced to the orthogonal one 
0(n,R) by the arguments \\xk — xl\\, K- [^Kt^Pl]- This may happen in such a 
way that V (2) is a sum of some purely affine term dependent only on T>, A4 and 
on an appropriate symmetry-restricting metrical term built of || • || and K,. 
It is a very interesting question whether the binary purely affine models 

V ( KL f = fKL(P,M) (79) 

may be realistic. The question was not yet touched seriously. Nevertheless, 
some limitations of applicability of the binary affine paradigm seem to be ob- 
vious. In our earlier papers O [10] we discussed dynamical models of a single 
affinely-rigid body, in particular, the purely geodetic models, i.e., ones without 
potentials. Lagrangian coincides then with the kinetic energies (metric ten- 
sors on GL(n,R)) given by (|TTj) or ([72)) , and first of all, by their special cases 
like ([73]) . ([74]). |75|) . It turns out that for incompressible affine bodies, when 
the configuration space of internal motion is restricted to SL(n,R), the purely 
geodetic affine models predict the existence of an open family of bounded (oscil- 
latory) trajectories within the general solution. However, on the non-restricted 
GL(n, R), when the volume changes are admitted, geodetic affine models predict 
the non-restricted dilatational motion, i.e., unlimited expansion or contraction. 
This is an evidently non-physical feature of these models. Therefore, at least 
some dilatations-stabilizing potential Vdii(det <p) must be assumed. When we 
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deal with systems of affine bodies, then it is clear that for an appropriate choice 
of / in (|79p the relative volumes 

det tp L / det ip K = det (ip K ~ 1 <PL,) = det T [ip K , ip L ] (80) 

are stabilized in the sense of performing bounded motions. However, no bi- 
nary potential may stabilize the single volumes det ipK themselves. Their time 
evolution will be non-bounded although the ratios (j80|) are bounded functions 
of time. To prevent this one should introduce some one-body potential term 
stabilizing (making bounded) the over-all dilatational behaviour, 

N 
K=l 

It is reasonable to assume that all V^J k are identical (when the body consists 
of identical elements) , 

N 

v£l(... ) <p A> ...)=J2f( det( PK)- (82) 

K=l 

(2) 

When V KL depend on their arguments in a proper way, so that &etT[ipK, <Pl] 
are bounded functions of time, then in principle it would be sufficient to use 
depending on det ip A for some fixed label A only. If such V^\detip A ) 
stabilizes det ip A , then automatically all volumes detipx will be stabilized by 
V^. But of course such a choice of the shape of V would not be either aesthetic 
or reasonable. 



3 General quantization ideas 

There is a direct logical chain from the atomic and molecular structure to macro- 
scopic properties, constitutive laws and material engineering. The point is par- 
ticularly delicate on the nano-level, where one is dealing with a very peculiar 
convolution of quantum and classical concepts. In any case, quantization is nec- 
essary then. Also some quasi-classical and correspondence problems are very 
relevant for these phenomena. 

The first step towards quantization is the classical canonical formalism [15[ 
[TT1 [8]. One should start from Legendre transformations which for potential 
systems with Lagrangians L = T — V(- • • ; Xk, <Pk] • • ■) are given by 

K _ dL dT Ka _ dL _dT_ 



or, alternatively, 

dv a K dv a K 



K dL dT ~ Ka _ ____ i v n 

P " 77v— L b - - - -3 (84) 



J 3 



dL 


dT 


dn K b a 


dn K " a 


dL 


dT 
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Inverting these formulas and substituting them to the energy expression 

FIT F)T 

E = v* K —+t K \j^--L (85) 

OV % K »6f a 

one obtains the classical Hamiltonian 

H = T + V. (86) 

Let us quote the resulting formulas for the geodetic (kinetic) Hamiltonians 
T. For the "usual" d'Alembert model © we obtain 

T d A = T* A + = _i_Tr (pp T ) + i T r (n T J~\) . (87) 

This is, as mentioned, the "usual" expression compatible with the d'Alem- 
bert principle. Although from some point of view it seems the best-motivated 
one, in complicated systems with collective modes and strong internal interac- 
tions some doubts and just objections may be raised against it. Our idea here 
was to concentrate on models motivated by symmetry principles, first of all, 
by afhne symmetry. Let us now review Legendre transforms of affine models 
quoted above. 

Of course, the model of translational kinetic energy T£ coincides 
exactly with Tf r A , so we have 

T tr is - af = —p T p=—TT(pp T ) 
2M P F 2M yFP ' 

= ^p T G[^p=^T r (pp T G[^), (88) 

just the first term of (|87p written in a few equivalent forms. 

The corresponding expression for (|65[) has the following form: 

T tl af - is = —p t p=—Tt(pp t ) 

Let us now quote the Legendre transforms of affinely-invariant internal ki- 
netic energies. For (|66|) we obtain 

= \^ (CM-V^J-V) = i T r (£ r J~i$) . (90) 

In particular, for the isotropic inertial tensor (|67|) leads to the expression: 

T^- is = ^Tr (£ t e) = -LTr (CM-V^tt) . (91) 

For the model (1681) we obtain 



3£"- = i T r (E^E) = i T r {fH tt) . (92) 



18 



Obviously, the Legendre transforms of the most general affine models (|7ip , (|72[) 
have the form: 



^7 af = ^C b a d c^ A b?. c d, (93) 



Vnt af = i^VfcS^E*,, (94) 

where the constant bimatrices £, 7?. are reciprocal to £, 7?. respectively. 

For the most interesting models (fT4")) . (|75|) , including their special doubly- 
affine case (|73l) wc obtain respectively 

^r af - ^ Tr ( sTs ) + ^ Tr ( s2 ) + ^( TrS ) 2 ' ( 95 ) 

zi 2^4 2i> 

1&+ = ^Tr(^ S ) + ^Tr(£ 2 ) + ^(Tr S ) 2 , (96) 

where the constants I, A, B are built of I, A, B in the following way: 

I = j(l 2 -A 2 ), A = j(A 2 -I 2 ), B = ~(I + A)(I + A + nB). (97) 



The special affine- affine case (f73|) corresponds to I = 0, and then obviously 
1/7 = and the first terms of |95|) . (|96|) do vanish. Of course, the second and 
third terms of (pj)) , (|9"6"|) are pairwise identical. 

These were expressions for various models of the kinetic Hamiltonians for 
single elements. Obviously, they should be labelled by the index K = 1, TV, and 
the total expression is obtained by a summation over K. 

Canonical formalism is very convenient and effective in analysis of classical 
equations of motion. They are represented then in terms of Poisson brackets, 

where F runs over some maximal functionally independent system of functions. 
To make use of (|98p one must establish the system of basic Poisson brackets 
for some geometrically distinguished quantities. As a rule, one uses quantities 
like S, E, because they are Hamiltonian generators of left and right regular 
translations in GL(n,K). Therefore, their Poisson brackets are determined by 
the structure constants of this group (and the mutual Poisson brackets between 
E and E-quantities do vanish of course, because the left translations commute 
with the right ones). Other important Poisson brackets are those between E, E- 
quantities and functions depending only on generalized coordinates. Calculating- 
such a bracket is identical with affecting configuration functions by first-order 
differential operators generating left and right regular translations in GL(n, K). 

However, the main advantage of Hamiltonian methods is that they provide 
a direct way towards quantization. 

Let us remind that the configuration space of our A^-body system is given 

by 

Q N * Q N tr x Q N iat ~ R nN x GL (n, M) N , 
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i.e., configurations are arrays ©. The manifold Q N is obviously an open subset 
of the linear space 

R nN x L (n, R) N ~ R nN x W l " N ~ R"("+ 1 ) JV . 

In any one-element configuration space R"x GL(n,R) we are given two distin- 
guished measures. One of them is the Haar measure a invariant under left and 
right group translations (cf. [9). The other one is the usual Lebesgue mea- 
sure a on R™x L(n,R). It is invariant under additive translations. In terms of 
coordinates 

da(x,(p) — dx 1 ■ ■ ■ dx n d(p 1 \ ■ ■ ■ d<p n n: (99) 
da (x, ip) — (det ip) n 1 da (x, ip) 

= (det ^r™" 1 da; 1 ■ ■ ■ dx n dip\ ■ ■ ■ d(p n n . (100) 

When we neglect translational motion then the Haar measure A on GL(n,R) 
and the Lebesgue measure I on L(n, R) are used 

dl(<p) = dip\---d<p n n , (101) 
d\(ip) = (detipy n ' 1 dl{ip) = (detipy n dip 1 1 ---dip n n . (102) 

Configuration spaces of the total iV-element system are endowed with the iV-told 
tensor products of these measures, a^ N \ a^ N \ 

The quantized theory is formulated in the following Hilbert spaces: 



L- 



(Q N ,aW) ,L 2 (g w ,aW) ,L 2 (GL(n,R) N A'")) , L 2 (GL(n, R) N , l^) 



Their elements, i.e., wave functions, are complex probability amplitudes of find- 
ing the system at a given classical configuration. Classical quantities depending 
only on configuration variables are represented in these L 2 -spaces as operators of 
multiplication by real- valued functions, in particular, by coordinates like x\ ip z a , 
etc. According to the general rules of quantum mechanics all other quantities 
are also represented by Hermitian or formally Hermitian (symmetric in dense 
domains) operators in these Hilbert spaces. Usually some ordering problems 
of non-commuting operators appear then. However, in dynamical applications, 
when Hamiltonian operators are constructed, one deals usually with some spe- 
cial physical quantities of well-defined geometric interpretation. As a rule, they 
are generators of symmetry groups underlying the problem. In our model they 
are just the affine spin in both the spatial and co-moving representation, the 
usual metrical spin and vorticity, etc. 

Linear momentum operators in spatial and co-moving representations are 
given respectively by 

Pa = 7^' Pq = 7^^' (103) 
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where, obviously, K — 1, TV is the "particle" label. These operators are formally 
Hermitian both in L 2 (Q N ,a^) and L 2 (Q N ,a ( - N '>). The operators 

X K \ = ?^ C -A-, % K * b = ^ K %—^ (104) 

l OipK c I OLpK a 

are formally Hermitian (not literally, they are unbounded as all differential op- 
erators) in L 2 (Q N ,a^) and in L 2 (GL(ra, R) N , A**') . Therefore, when using 
these Hilbert spaces we may interpret (I104|) as operators of affine spin respec- 
tively in the spatial and co-moving representations. 

Just as in classical theory, p K a are infinitesimal generators of translations of 
the K-th constituent. Similarly, T,K a b generate spatial affine transformations 
(rotations and homogeneous deformations) of internal degrees of freedom of the 
K-th "molecule". T,K a b generate material affine transformations of the if-th 
element. Namely, let us consider the operators 

V K (y) := expQyV^a), V G R", (105) 

L K (z) := exp^z b a S K a ^, z G L (n, M.) , (106) 

where the operator exponent is meant in the usual power-series sense. If this 
series convergent in the action on some function 'J : Q N — > C, then 

{V K (y)9)(...,XA,...;...,<p B ,...)= (107) 
= ,x A + yS AK , ip B , ...) , 

(L K (z)^)(...,x A ,...;...,^ B ,...) = (108) 
= *(..., x A , exp (zJkb) <^b, • • •) • 

Similar statements may be formulated about the co-moving objects, e.g., defin- 
ing 

(z) := exp (^z b a 'Z K a b ^ , z G L (n, R) , (109) 

we obtain that 

(R K (z)y)(...,x A> ...;...,tp B ,...) = (110) 
= (..., xa, tfB exp (zSbk) ,■■■)• 

One can act separately on all arguments, nevertheless, the special geometric 
role is played by transformations acting in the same way on all arguments, e.g., 

V( y ) = Vi(y)---V N (y), (Hi) 

L(z) = U(z)---Ln(z), (112) 

r(z) = n 1 (z)---n N (z). (ii3) 



21 



Their generators are respectively identical with the total linear momentum and 
the total afhnc spin in the spatial and co-moving representations, 



N N N 



Pa = J2 p K -> SQ& = S ^ a "' ga& = ^ a «- ( 114 ) 



K=l K=l K=l 

Obviously, 



V(y) = exp(-y a Pa ), (115) 



L(z) = exp^z b Q S a h J, (116) 

R(z) = expQz\£ a ^. (117) 

Obviously, all exponential operators quoted here are unitary in L 2 (Q N ,a^ N ') 
or L 2 (GL(n,R) Ar , A (JV ^). However, L K , R K are not unitary in L 2 (Q N ,a^) 
and L 2 (GL(n,R) JV ,^ Ar )). The reason is that the measures a, I are not in- 
variant under group translations. Also, when working in L 2 (Q N , and 
L 2 (GL(n, M) N , l( N A , that is admissible, one must modify the definition of the 
above unitary operators (introducing some multipliers). The generators (|104l) 
are not formally Hermitian and to become such they must be modified by some 
additive corrections: 

Tun --s ^-v. Tun 

^K a b '■= ^K a b + 7TT<5 a f), Sif °b := Sjf"(, + — S a b- (118) 
ll ll 

Let us also quote the formally Hermitian operators 

J K \ = XK a p K b + V K a b = + V 6 , (119) 

which generate affine transformations acting both on translational and internal 
degrees of freedom of the K-th constituents. Ak and Y,k are respectively the 
translational (orbital) and internal parts. One can also introduce the total 
quantities 

J a 6 = A a 6 + S Q b (120) 

obtained by the if -summation. 

In analogy to (|118p . we have that 

'A K a b := K K \ + ^5\. (121) 
li 

Let us observe that for the total quantities built of (|1 18[) we have 

= E« 6 + ^S a bl '£% = E« 6 + ^S\. (122) 
li li 
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and similarly 

'A\ = A\ + —6\. (123) 
2^ 

After quantization the canonical momenta ir Ka i conjugate to tfK l a become 
operators: 

k-„ h d , 

They are formally Hermitian in L 2 (Q N ,a^), L 2 (GL(n, M.) N , jW), but not in 
L 2 (Q'f.aW), L 2 (GL(n,R) Ar , A (Ar) ), so now the situation is quite opposite to 
the previous one. 

Hamilton operator has the following form: 

H = T + V, (125) 

where T is the kinetic energy operator and V is the potential term; usually one 
does not distinguish graphically between the function V and the operator V 
which multiplies the wave function by V, 

(y*)(x,<p):=V(x,<p)9(x,<p). (126) 

In usual structural problems quantum dynamics reduces to the stationary Schro- 
dinger equation, i.e., to the energy following eigcnproblcm: 

= Ety. (127) 

And finally, we see that constructing the kinetic energy operator T is a crucial 
step of the quantization procedure. And this may be done just on the basis of 
the above classical expressions for kinetic energies in terms of canonical phase- 
space variables ftT), ftg), (HJ, (JgDJ), (EE]), (M}, (S3), (ED, (G3, HMD- Simply 
one should algebraically substitute the above operators p K i, p K a , P Ka j, ^ Ka b, 
Y, Ka b instead of the corresponding classical expressions. Because of the geomet- 
ric meaning of these quantities as generators of natural transformation groups, 
there is no problem of ordering of operators. The only point one should be care- 
ful with is just the one concerning the measures used in configuration spaces 
when defining the L 2 -Hilbert spaces. In affine models the measures a, A are 
more natural and then one uses the purely differential operators (|103[) , (|104|) . 
In Hilbert spaces based on a, I we would have to use the modified expressions 
(|118|) . And conversely, in non-affine d'Alembert models the Lebesgue measures 
a, I are more natural, because they enable one to use the purely differential 
operators (] 1 24[) without any algebraic correction. 

The above geometric approach is very convenient. If we tried to calculate 
the kinetic energy operators as 

= 2 

where A is the d'Alembert operator based on the metric tensor underlying 
the classical expression T, the calculations would be hopeless and the result 
completely obscure, non-useful. 
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When solving any particular problem one must use coordinates. For our 
purposes the most convenient choice is that based on the polar and two-polar 
decompositions: 

ip = UA = BU = LDR- 1 = LDR T , 

where <p £ GL+(n,R), U,L,R e SO(n,M) (orthogonal), A, B = UAIJ- 1 are 
symmetric and positively definite, and D is diagonal and positive. 

Green and Cauchy deformation tensors are then expressed as follows: 

G = A 2 = RD 2 R T , C = B- 2 =LD- 2 L-\ 

It is convenient to denote 

D aa = Q a = cxp (<?«). 

The quantities Q a , q a offer another convenient choices of basic deformation 
invariants. The Haar and Lebesgue measure A, I are then expressed as 

d\(L,D,R) = n|sh(V -q^)\dq 1 ...dq n dfi(L)d t i{R), 

di (l, d,r) = Y[ (Q l - Q°) (Q l + Q j ) d Q'--- d Q nd v ( L ) d v ( R ) . 

where \x is the Haar measure on SO(n,R). Performing some partial integra- 
tions one can obtain from probability distributions for some quantities like 
deformation invariants, orientations of the main axes of Green and Cauchy defor- 
mation tensors and also probability distributions for particular values of Green 
and Cauchy deformation tensors. 
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